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We present a calculation of the nucleon sigma term on two-flavor QCD configurations with dy- 
namical overlap fermions. We analyse the lattice data for the nucleon mass using the baryon 
chiral perturbation theory. Using partially quenched data sets, we extract the connected and dis- 
connected contributions to the nucleon sigma term separately. Chiral symmetry on the lattice 
simplifies the determination of the disconnected contribution. We find that the strange quark con- 
tent, which determines the neutralino dark matter reaction rate with nucleon through the Higgs 
boson exchange, is much smaller than the previous lattice results. 
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1. Introduction 

Nucleon sigma term Gjin is given by a scalar form factor of nucleon at zero recoil. While up 
and down quarks contribute to Gjin both as valence and sea quarks, strange quark appears only as a 
sea quark contribution. As a measure of the strange quark content of the nucleon, the y parameter 
is commonly introduced. These parameters are defined as 

2(N\ss\N) 

OnN = mud{N\uu + dd\N), y= \ ' ' . (1.1) 

{N\uu + dd\N) 

The y parameter plays an important role to determine the detection rate of possible neutralino dark 
matter in the supersymmetric extension of the Standard Model [|l|, ||, ^ ||]. Already with the 
present direct dark matter search experiments one may probe a part of the MSSM model param- 
eter space, and new experiments such as XMASS and SuperCDMS will be able to improve the 
sensitivity by 2-3 orders of magnitude. Therefore, a precise calculation of the y parameter will be 
important for excluding or proving the neutralino dark matter scenario. 

Using lattice QCD, one can calculate the nucleon sigma term directly. Furthermore, it is 
possible to determine the valence and sea quark contributions separately. Previous lattice results 
were o^n = 40-60 MeV, y = 0.66(15) [|], and o^n = 50(3) MeV, y = 0.36(3) [0] within the quenched 
approximation, while a two-flavor QCD calculation ^ gave Ojin = 18(5) MeV and y = 0.59(13). 
There is an apparent puzzle in these results: the strange quark content is unnaturally large compared 
to the up and down contributions that contain the connected diagrams too. 

Concerning this problem, it was pointed out that using the Wilson-type fermions, the sea 
quark mass dependence of the additive mass renormalization and lattice spacing can give rise to 
a significant lattice artifacts in the sea quark content [^. Unfortunately, after subtracting this 
contamination the unquenched result has large statistical error, y = —0.28(33). 

In this study, we analyze the data of the nucleon mass obtained from a two-flavor QCD simula- 
tion employing the overlap fermion [ [lO| ] which can remove this problem by explicitly maintaining 
exact chiral symmetry on the lattice. Although the two-flavor QCD calculation cannot avoid the 
systematic error due to the neglected strange sea quarks, our study with exact chiral symmetry re- 
veals the underlying systematic effects in the calculation of the nucleon sigma term, especially in 
the extraction of its sea quark contribution. It therefore provides a realistic test case, which will be 
followed by 2-1-1 -flavor calculations in the near future The full details of this work is presented 



2. Lattice simulation 

We make an analysis of the nucleon mass obtained by two-flavor QCD configurations gen- 
erated with dynamical overlap fermions. Our simulations are performed at a lattice spacing a = 
0.118(2) fm on a 16^ x 32 lattice with a trivial topological sector 2 = 0. For each sea quark mass, 
we accumulate 10,000 trajectories; the calculation of the nucleon mass is done at every 20 trajec- 
tories, thus we have 500 samples for each mgea- For the sea quark mass amsea vve take six values: 

'For a very recent result from 2+1 -flavor QCD, see |p^. 
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0.015, 0.025, 0.035, 0.050, 0.070, and 0.100 that cover the mass range mj6-ms with m^. the phys- 
ical strange quark mass. Analysis of the pion mass and decay constant on this data set is found in 



[]13[]. In order to improve the statistical accuracy, we use the low-mode preconditioning technique. 
The two-point function made of low-lying modes of the overlap-Dirac operator is averaged over 
different time slices with 50 chiral pairs of the low modes. For the quark propagator, we take a 
smeared source defined by a function 0(|^|) oc exp(— A|3c|) with a fixed A = 0.40. We then calculate 
the smeared-local two-point correlator and fit the data with a single exponential function after av- 
eraging over forward and backward propagating states in time, we take the valence quark masses 
amvai = 0.015, 0.025, 0.035, 0.050, 0.060, 0.070, 0.080, 0.090, and 0.100. 

The matrix element defining the nucleon sigma term can be related to the quark mass de- 
pendence of the nucleon mass using the Feynman-Hellman theorem, which derives the relations 
between the nucleon mass and the quark contents of the nucleon as 

^^ = {N\{uU + dd)\N),onn. ={N\{uU + dd)\N) disc- (2.1) 

The subscripts "conn" and "disc" on the expectation values indicate that only the connected or 
disconnected quark Une contractions are evaluated, respectively. In the present study we exploit 
this indirect method to extract the matrix elements corresponding to the nucleon sigma term. 

Another possible method is to directly calculate the matrix element from three-point functions 
with an insertion of the scalar density operator {Uu + dd){x). In principle, it gives a mathematically 
equivalent result to the indirect method including lattice artifacts. Numerical differences could 
arise only in the statistical error and the systematic uncertainties of the fit ansatz. 

3. Analysis of the unitary points with baryon chiral perturbation theory 



We carry out baryon chiral perturbation theory (BChPT) [14] fits of the nucleon mass at the 



five heaviest quark masses for unitary point by using a simplified fit function 



Ma, = Mo - 4ciml - ^l^'"'^ + e[ {l^)ml, (3. 1) 

where Mq is the nucleon mass in the chiral limit and fji is the pion decay constant fixed at its 
physical value 92.4 MeV and the constant gA describes the nucleon axial-vector coupling. Since 
the coupling is very well known experimentally, we attempt two options: (Fit a) a fit with fixed 
gA (=1.267), and (Fit b) a fit with gA being dealt as a free parameter. Moreover we also carry out 
the full ^{p^) fit and ^{p^) fit following the analysis done in [15]. The left panel of Fig. |l| shows 



the BChPT fits. The lattice data show a significant curvature towards the chiral limit. All the details 



of our analysis can be found in []1 IQ. 

In order to estimate the systematic uncertainties due to the finite volume effect , we correct the 
data for the finite volume effect using BChPT &{p'^) formula [16]. For the input parameters Mq, 



gA, c\, we use the nominal values (Mq = 0.87 GeV, gA = 1.267, ci = —1.0 GeV ). The chiral fit 



is then made for the corrected data points using the formula (3.1) for the five or six heaviest data 
points. The result is shown in the right panel of Fig. After correcting the finite volume effect, 
there are 5-8% decrease in Mq and 4-7% increase in the magnitude of the slope |ci|. We observe 
that the deviation due to the finite volume effect is smaller than the uncertainty of the fit forms. 



3 



Nucleon sigma term and strange quark content with two flavors of dynamical overlap fermions H. Ohki 




m^' [GeV"] m^' [GeV"] 



Figure 1: Chiral fit of the lattice raw data (left) and the corrected data (right). In the left panel, the solid 
curve represents the Fit a. The dot, dashed, dot-dashed curves correspond to the fits using BChPT including 
higher order terms. In the right panel. Solid and dashed curves represent the fits using 5 and 6 heaviest data 
points (diamond), respectively. For a reference, we also show the raw data (circles). 

Taking the Fit a (gA fixed, Finite volume corrections not included) as our best fit and using the 
variation with fit ansatz and Finite volume corrections as an estimate of the systematic errors, our 
result of the nucleon sigma term is 

a^N = 52(2),tat(^'7 )extrap(^o)FVE McV, (3.2) 

where the errors are the statistical and the systematic due to the chiral extrapolation (extrap) and 
finite volume effect (FVE). The largest uncertainty comes from the chiral extrapolation. The fi- 
nite volume effect is sub-leading, which is about 9%. This result is in good agreement with the 
phenomenological analysis based on the experimental data. 

4. Partially quenched analysis and extraction of y 

In order to extract the connected and disconnected diagram contributions separately, we fit the 
quark mass dependence of the nucleon mass with the one-loop partially quenched ChPT formula 


Mn = Boo + Bio(<f +BoiK^)' +B2o(<)' +Bn(<f +Bo2(<)' 
'l6^{i [-n<? + m<?+9m-\m^^f)]+^[-l9{nr^;)'^ 

+ ^ [_13«)3 +4(,<y +9m-«)')] } (4.1) 

where m''^', m™, and m^* denote the pion mass made of valence- valence, valence-sea, and sea-sea 
quark combinations, respectively. At this order, one can rewrite the expression in terms of mvai and 
wzsea using the leading-order relations. The coupling constant gi is another axial-vector coupling. 
We use nominal values gA = 1.267 and g\ = —0.66 when they are fixed in the fit. There are also 
contributions from the decuplet baryons. In our analysis these contributions can be absorbed into 
the analytic terms in (pTT|). The independent fit parameters are Boo> ^oi> 5io> fiii> B20, Bo2> gi, and 
gA- We attempt a fit with fixed gA and gi (Fit PQ-a), a fit with fixed gA (Fit PQ-b) and a fit with all 
the free parameters (Fit PQ-c). Finite volume corrections are not taken into account. 
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Figure 2: The left panel shows partially quenched nucleon masses and fit curves (Fit PQ-b). The right panel 
shows the connected (up) and disconnected (down) contributions to the sigma term evaluated at myai — m^^a.- 
Solid, dashed and dotted curves represent the results from the Fit PQ-a, PQ-b and PQ-c, respectively. The 
two vertical lines show physical up and down quark mass aniq ~ 0.0034 (left) and strange quark mass 
am^ = 0.084 (right). 

The left panel of Fig. ^ demonstrates the result of the partially quenched ChPT fit. The sea 
quark mass dependence at eight fixed valence quark masses are nicely fitted with the formula (p^). 
We find that both axial-couplings can be determined with reasonable accuracy {gA = 0.93(22) and 
gi = —0.29(5) ) for the Fit PQ-c. The right panel of Fig. § shows the partial derivatives in ( ^Tl] ) 
with respect to myai and to mgea evaluated at the unitary points myai = /Wsea- For both contributions, 
we clearly find an enhancement towards the chiral limit. Results with different fit ansatz show 
slight disagreement near the chiral limit, which indicate the size of the systematic uncertainty. We 
find that the sea quark content of the nucleon {N\{uu + dd)\N)iif,i^/ {N\{uu + dd)\N)i^ona is less than 
0.4 for the entire quark mass region in our study, so that the valence quark content is the dominant 
contribution to the sigma term. 

Rigorously speaking, it is not possible to extract the strange quark content {N\ss\N) within 
two-flavor QCD. Instead, in this work, we provide a "semi-quenched" estimate of the strange 
quark content which is defined as the ratio of the strange quark content (discomiected contribu- 
tion at nival = "^sea = wJ.v) to the up and down quark contributions (connected plus disconnected 
contributions at myai = ?Msea = mud) following [^. Taking the result from the Fit PQ-b as a best 
estimate, we obtain the parameter y as 

//=2 = 0.030(16)stat(i8)extrap(t2)"..' (4-2) 

where the errors are statistical and systematic from chiral extrapolation and from the uncertainty 
of rus, respectively. The chiral extrapolation error is estimated by the differences of the results of 
Fit PQ-a,b and c. We also note that there may be an additional ~10% error from the finite volume 
effect as discussed in Section ||, but it is much smaller than the statistical error in our calculation. 



5. Discussion and summary 

We found that the disconnected contribution to the sigma term is much smaller than the previ- 
ous lattice calculations with the Wilson-type fermions y ~ 0.36 ~ 0.66 [|[, ^, (except for as 
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Figure 3: Two-loop diagrams for the mass shift and the operator mixing involving both the valence and 
the sea quarks. These graphs show the corresponding contribution to the sigma term. Due to the power 
divergent mass shift, the flavor singlet sea quark scalar operator insertion contributes as the flavor singlet 
scalar operator insertion with coefficient of additive mass shift. 



explained below). The authors of found that the naive calculation with the Wilson-type fermions 
may over-estimate the sea quark content due to the mixing effect from the additive mass shift as 



^M^ 



„bare 
"val 



^M^ 



+ 



1 phys 
dm. 



'val 



„phys 
"val 



^^bare 



„bai-e 
'^val 



(5.1) 



where Z,„ is the renormalization factor of the quark mass. Therefore, in order to obtain the deriva- 



tive (gjj) one must subtract the unphysical contribution from the additive mass. They found that 
their unsubtracted result y = 0.53(12) is substantially reduced and becomes consistent with zero: 
y = —0.28(33). It should be noted that this mixing effect is not limited to the spectrum method 
but arises also in the direct matrix element method. It is easy to see that the operator is indeed a 
derivative of the mass shift with respect to the sea quark mass so that 



bare 
sea 



(w) 



phys 



+ 



3 phys 
am. 



'val 



/val 



,,bare 
"val 



(5.2) 



holds. (Subtraction of the divergent counter terms is assumed on the left hand side.) Fig. || shows 
the two-loop contribution of the operator mixing of the bare sea quark operator with physical va- 
lence quark operator due to the additive mass shift which arises from the absence of the chiral 
symmetry. Our calculation using the overlap fermion is free from this artifact because of the exact 



chiral symmetry. Therefore, the small value of y obtained in our analysis ( ft.2D provides a much 
more reliable estimate than the previous lattice calculations ^. 

In summary, we have calculated the nucleon sigma term and the strange quark content of 
nucleon in two-flavor QCD simulation on the lattice with exact chiral symmetry. Owing to the 
exact chiral symmetry, our lattice calculation is free from the large lattice artifacts coming from the 
additive mass shift present in the Wilson-type fermion formulations. From an analysis of partially 
quenched lattice data, we have found that the sea quark content of the nucleon is less than 0.4 for 
the entire quark mass region in our study. 

An obvious extension of this work is the calculation including the strange quark loop in the 
vacuum. Simulations with two light and one strange dynamical overlap quarks are on-going [[l9|]. 



^Strictly speaking, there is another artifact from the sea quark mass dependence of the lattice spacing. We have found 
that the lattice spacing dependence is negligibly small, since we exploit mass-independent renormalization scheme (c.f. 
mM) in contrast to \W, in which they used the Sommer scale for the scale input at each sea quark mass. 
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The main numerical calculations were performed on IBM System Blue Gene Solution at High 
Energy Accelerator Organization (KEK) under support of its Large Scale Simulation Program 
(No. 07-16). We also used NEC SX-8 at Yukawa Institute for Theoretical Physics (YITP), Kyoto 
University and at Research Center for Nuclear Physics (RCNP), Osaka University. The simulation 
also owes to a gigabit network S1NET3 supported by National Institute of Informatics for efficient 
data transfer through Japan Lattice Data Grid (JLDG). This work is supported in part by the Grant- 
in- Aid of the Ministry of Education (Nos. 18340075, 18740167, 19540286, 19740121, 19740160, 
20025010, 20039005, 20740156). The work of HP is supported by Nishina Memorial Foundation. 
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